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A new mechanism for the enhanced generation of compressible fluctuations by Alfven waves is 
presented. A strongly nonlinear regime of Alfven wave phase-mixing is numerically simulated in a 
one-dimensionally inhomogeneous plasma of finite temperature. It is found that the inhomogeneity 
of the medium determines the efficiency of nonlinear excitation of magnetoacoustic waves. The level 
of the compressible fluctuations is found to be higher (up to the factor of two) in inhomogeneous 
regions. The amplitude of the generated magnetoacoustic wave can reach up to 30% of the source 
Alfven wave amplitude, and this value is practically independent of the Alfven wave amplitude 
and the steepness of Alfven speed profile. The highest amplitudes of compressible disturbances 
are reached in plasmas with /3 of about 0.5. The further growth of the amplitude of compressible 
fiuctuations is depressed by saturation. 
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Magnetohydrodynamic (MHD) waves play an impor- 
tant role in the dynamics of laboratory, space and astro- 
physical plasmas (e.g. [Q). The simultaneous existence 
of several MHD wave modes, three, at least, in the sim- 
plest case, stimulates the interest to processes of mode 
coupling and interaction. In particular, the interaction 
of linearly incompressible Alfven waves with compress- 
ible magnetoacoustic waves is essential in problems such 
as the heating of the open corona of the Sun and accel- 
eration of the solar wind. Indeed, Alfven waves are often 
named as the carriers of the energy from the lower layers 
of the solar atmosphere to the corona. Also, in the inner 
heliosphere, Alfven waves represent the main component 
in MHD turbulence 1^ . However, compressible waves are 
subject to more efficient dissipation, as they decay on 
volume viscosity, not shear as to incompressible Alfven 
waves, and the difference of these two viscosities can be 
several orders of magnitude. Moreover, this mechanism 
allows energy to be transported across field lines in con- 
trast to Alfven waves which can only transport energy 
along the field. Also, compressible waves perturb the 
plasma density and consequently can be detected, again 
in contrast to Alfven waves, with imaging telescopes. 

A classical example of such interaction is the decay in- 
stability of Alfven waves, connected with the resonant 
three-wave interaction of Alfven and magnetoacoustic 
waves 1^. The efficiency of this phenomenon is deter- 
mined by the amplitudes of the interacting waves. How- 
ever, this mechanism works only for quasi-periodic (per- 
haps wide-spectrum, [Q) waves, and is not so prominent 
for wave pulses that could be generated e.g. by some ex- 
plosive events such as solar flares, coronal mass ejections, 
etc. 

In contrast, non-resonant mechanisms for compress- 
ible fluctuation excitation are independent of the wave 
spectrum and can be efficient even for short wave pulses. 
An example of non-resonant MHD wave interaction is 
the nonlinear excitation of longitudinal magnetoacoustic 



perturbations by nonlinear elliptically polarized Alfven 
waves through the ponderomotive force. This phe- 
nomenon is actually the generation of the second har- 
monic and, in contrast with three wave resonant inter- 
action, and does not require the daughter wave to be 
present in the system from the very beginning. Also, 
the efficiency of this process is restricted by the source 
wave amplitude, but does not require the source wave 
to be harmonic. Another parameter affecting the ef- 
ficiency is /3, the thermal to magnetic pressure ratio. 
The generation of the compressible perturbations leads to 
self-interaction and consequent steepening of the Alfven 
wave, whose evolution is described by the Cohen-Kulsrud 
equation |5|. 

Yet another non-resonant mechanism for the genera- 
tion of compressible fluctuations arises when the Alfven 
speed varies across the magnetic field. Initially plane, lin- 
early polarized Alfven waves become oblique and sharp 
gradients in the direction across the field are secularly 
generated. This phenomenon is known as Alfven wave 
phase mixing, and it has been intensively studied in the 
context of the solar coronal heating problem l6| and of- 
ten seen in full MHD numerical simulations (e.g. [Q, 
Q). In the compressible regime, phase mixing of Alfven 
waves is accompanied by nonlinear generation of fast 
magnetoacoustic waves In this phenomenon, the ef- 
ficiency of nonlinear coupling is dramatically affected by 
the inhomogeneity of the medium. To illustrate this, con- 
sider a low-/? plasma with the straight and uniform mag- 
netic field Bq = BqBz, where Bq is the absolute value 
of the unperturbed field and is the unit vector in z- 
direction. The plasma mass density pq{x) varies across 
the field, in the x-direction. Thus, the Alfven speed 
CAix) — Bq[4t:pq{x)]^^^^ also varies across the field. In 
this geometry, Alfven waves (perturbing Vy and By) and 
fast magnetoacoustic waves (perturb Vx, B^, B^ compo- 
nents of the bulk velocity, magnetic field and the mass 
density p) are linearly coupled if the Alfven waves have 
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a component of k in the y-direction [Q. However, when 
d/dy tends to zero, the waves become linearly decoupled 
from each other and in this case, the effect of nonhn- 
ear coupling becomes pronounced. The initial stage of 
fast magnetoacoustic wave generation by transverse gra- 
dients in weakly nonlinear Alfven waves is described by 
the equation 



-B, 



dt^ ^ ^ ^ ' \ dx^ ~^ dz^ 87TP0 (x) dxdz ■ 

_ (1) 

The right handside of Eq. in the Fourier transform 
domain is proportional to k^kxA^, where kz and k^ are 
the components of the wave number in the longitudi- 
nal and transverse directions respectively and A is the 
Alfven wave amplitude. Consequently, even in the case 
of homogeneous Alfven speed {Ca — const) ^ a non-plane 
{kx ^ 0) Alfven wave generates compressible perturba- 
tions, which, according to the left hand side of Eq. (|^), 
propagate isotropically. 

In the weakly nonlinear regime the Alfven wave steep- 
ening can be neglected and the linear solution can be 
used, 
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FIG. 1: Contourplots of the evolution of an initially plane 
single Alfvenic pulse in a plasma with an Alfven speed varing 
in the a;- direct ion. The unperturbed magnetic field is straight 
and has the z-component only. Here A = 0.5, A = 0.31, /? = 2. 



By = Af[z-CA{x)t\, 



(2) 



where / is an arbitrary smooth function prescribed by 
the initial conditions. 

With an inhomogeneous Alfven speed, the Alfven wave 
is subject to phase mixing and kx ^ oo as time pro- 
gresses. Consequently, the phase mixing amplifies the 
right hand side term of Eq. (|^) and affect MHD wave 
coupling. Indeed, with the Alfven wave given by Eq. (||), 
the right hand side of Eq. becomes 
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and, consequently, it grows secularly in time, proportion- 
ally to the product of the Alfven wave amplitude and 
the inverse characteristic spatial scale of the medium in- 
homogcneity (see |9j, for details). The secularity in 
the nonlinear generation of the fast waves is entirely con- 
nected with the inhomogeneity of medium. Thus, even if 
the Alfven wave amplitude is weak, the generation of fast 
wave can be dramatically amplified by the medium inho- 
mogeneity. In the homogeneous case, the same pondero- 
motive force can also generate perpendicular compress- 
ible fluctuations if the Alfven wave is initially non-plane. 
However, in this case, the perpendicular gradients in the 
Alfven wave remain as prescribed by the initial condi- 
tions and the generation of compressible fluctuations is 
limited by the initial shape of the wavefront. The induced 
longitudinal motions, also generated by the Alfven wave, 
and do not grow secularly either. 

The secular generation of fast magnetoacoustic waves 
by phase-mixed Alfven waves in initial stage of the wave 
interaction was observed in fuU-MHD numerical simula- 
tions [p|, ItI]. The adequate modelling of the developed 




FIG. 2: A snapshot of the density perturbation at t = 15, 
generated by an Alfvenic pulse of the relative amplitude 0.5, in 
a plasma with (3 = 2 and inhomogeneity parameter A = 0.31. 



stage of this phenomenon could be performed with very 
high numerical resolution that became possible only re- 
cently. Indeed, it was found ||ll| that in the weakly non- 
linear regime, phase-mixing leads to quick saturation of 
the interaction. In the regime, the saturation does not 
allow the fast wave amplitude to grow to more than a 
few tenth of a percent of the initial Alfven wave am- 
plitude. However, the tendency of the saturation level 
to grow with the Alfven wave amplitude was noticed, 
which indicated that in the case of higher amplitudes, 
the generation of compressible waves can be significant. 
In the strongly nonlinear regime, in addition to phase 
mixing, coupling of MHD waves is also affected by the 
steepening of the Alfven wave j5| , and by the back reac- 
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FIG. 3: Evolution of maximal values of the relative ampli- 
tude of the generated compressible perturbations, transverse 
max{\Vx{x, z,t)\) and longitudinal maa;(|14(x, 2, and that 
of Alfven wave 'max{\Vy{x,z,t)\), in time. The solid curve 
presents decay of the initial Alfven perturbation due to shock 
dissipation in the case of no plasma density inhomogene- 
ity (A — 0). The short dashed curve with triangles de- 
picts max{\Vy{0, z,t)\) when A — 0.31, while dotted curve 
with astersics represents the same physical quantity on the 
edge of the similation box (away from the phase-mixing re- 
gion) i.e. max {\Vy {15, z,t)\). The dash-dotted curve with 
diamonds represents max{\Vz{x, z,t)\), while the dash-triple- 
dotted curve shows max{\Vx{x, z,t)\) both for the case of 
A = 0.31. The long dashed curve represents max{\Vz{x, z, t)\) 
for the case when A = 0. Here A = 0.5 and (3 — 2. 

tion of the generated fast wave on the Alfven wave. In- 
deed, as it has been shown by numerical simulations 0, 
the nonlinear interaction of magnetoacoustic and Alfven 
waves causes energy transfer to smaller and smaller spa- 
tial scales. This phenomenon may lead to the forma- 
tion of shock waves and, therefore, to increased heating 
of plasma by viscosity and resistivity. In modelling of 
non-resonant interaction of MHD waves on a plasma in- 
homogeneity, we consider a plasma configuration similar 
to the one investigated in 0| and i.e. the plasma has 
one-dimensional inhomogeneities in the equilibrium den- 
sity po{x) and temperature To{x), and is penetrated by a 
straight and homogeneous magnetic field directed along 
the z-axis. The unperturbed total pressure is taken to 
be constant everywhere. The density profile is a smooth 
interface, 

3 - 2tanh(Aa;), (4) 

where A is a parameter prescribing the steepness of the 
profile near a; = 0. The density is normalized to po{x — 
0). The sharpest gradients are located near x — 0. Both 
the Alfven speed CAix) and the sound speed Cs{x) are 
inhomogeneous in the x-direction, however the param- 
eter (3 and the unperturbed total pressure are constant 
across the profile. Dynamics of the plasma is described 
by single-fiuid MHD equations in Cartesian coordinates. 



FIG. 4: Top panel: Maximal values of the relative ampli- 
tude of the generated compressible perturbations for different 
values of the inhomogeneity parameter A, with /3 = 2 (solid 
curve) and l3 = 0.5 (dashed curve). The initial amplitude of 
the Alfvenic pulse is ^4 = 0.5. Mid panel: Maximal values of 
the relative amplitude of the generated compressible pertur- 
bations for different values of initial amplitude of the Alfvenic 
pulse, with inhomogeneity parameter A = 0.5 and (3 = 1.25. 
Bottom panel: Maximal values of the relative amplitude of 
the generated compressible perturbations for different values 
of /3, with initial amplitude of the Alfvenic pulse A = 0.5 and 
inhomogeneity parameter A = 0.5. 



The MHD equations are solved with Lare2d |T^. 

An initially plane linearly polarized Alfvenic pulse of 
Gaussian shape and highly nonlinear amplitude Vy = 
0.5Ca{0), propagates along the background magnetic 
field, in the z-direction. There are no compressible per- 
turbations present in the system at the beginning. The 
pulse is subject to phase-mixing, which takes place in 
the vicinity of a; = 0. Fig. 1 shows the development 
of the transverse gradients in the Alfven pulse. Longi- 
tudinal and transverse gradients generate compressible 
perturbations {Vx, Vz, p, and B^), associated with 
slow and fast magnetoacoustic waves. Note, that initially 
plane Alfven waves (perturbing Vy and By) in the homo- 
geneous medium do not generate the perturbations of Vx , 
Bx and Bz. Thus, the generation of these perturbations 
in our simulations are caused by the effect of plasma in- 
homogeneity (see Fig. 3). It is clearly seen (Fig. 2) that 
the relative perturbations of the plasma density are en- 
hanced by a factor of about two in the phase-mixing re- 
gion (near a; = 0). As in the weakly-nonlinear regime 
( [prl ), initially the compressible disturbances grow secu- 
larly, according to Eq. (^), and then a saturated state is 
approached (Fig. 3). The saturation time is proportional 
to the efficiency of the generation. However, the maxi- 
mal relative amplitude of compressible perturbations is 
practically independent of the problem paramteres. The 
parametric study (Fig. 4, top and mid panel) demon- 
strates that the efficiency of the non-resonant genera- 
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tion of compressible disturbances by phase-mixed Alfven 
waves depends weakly upon the initial Alfven wave am- 
plitude and the steepness of the inhomogeneity, repre- 
sented by the parameter A (see Eq. (^)). The depen- 
dence on the plasma parameter (3 is more pronounced 
(see Fig. 4, bottom panel), and the highest amplitude of 
compressible disturbances is reached for (3 ~ 0.5. For 
all investigated parameters, the observed saturation level 
was about 30% of the initial Alfven wave amplitude. 

We conclude that in the presence of plasma inhomo- 
geneity across the magnetic field, incompressible (Alfven) 
and compressible (magnetoacoustic) MHD fluctuations 
are effectively coupled by non-resonant mechanisms. 
This coupling takes place even if the waves are linearly 
decoupled, e.g. the perturbations are plane in the ignor- 
able direction perpendicular to both the magnetic field 
and the inhomogeneity gradient. These results demon- 
strate that the efficiency of the nonlinear interaction of 
the different type modes is strongly affected by the inho- 
mogeneity of the medium. In particular, in the presence 
of a one-dimensional inhomogeneity, nonlinear generation 
of compressible disturbances by incompressible Alfven 
waves is enhanced up to a factor of two (Fig. 3). The am- 



plitude of the generated compressible disturbances expe- 
riences saturation at the level of about 30% of the initial 
amplitude of the source Alfven wave. For example, it is 
10-20% of the background density for the Alfven wave 
amplitude of about 0.5. In particular, this result is rel- 
evant to the explanation of the absence of a significant 
compressible component in the solar wind MHD turbu- 
lence. Indeed, the density fluctuations of less than 30% 
of the background value, prescribed by the saturation, 
can hardly be detected in a high-/3 plasma because of a 
high level of the thermal noise in the in situ data. These 
results are especially relevant to the interpretation of ob- 
servational data obtained by CLUSTER-II mission. This 
issue will be discussed in more detail elsewhere. 

Finally, we would like to point out that the pro- 
posed mode coupling mechanism could find applications 
in other branches such as propagation of elastic waves 
in a inhomogeneous, fluid-saturated porous medium, 
as mathematically, the equations describing this phe- 
nomenon are similar to the MHD equations p^ . 

Numerical calculations of this work were done using 
the PPARC funded Compaq MHD Cluster in St An- 
drews. 



[1] M. Goossens, Space Sci. Rev. 68, 51 (1994); B. Roberts, 

Solar Phys. 193, 139 (2000). 
[2] C.-Y. Tu and E. Marsch, Space Sci. Rev. 73, 1 (1995). 

B.T Tsurutani and CM. Ho, Rev. Geophys. 37, 517 

(1999). 

[3] A.A. Galeev and V.N. Oraevskii, Sov. Phys. Dokl. 7, 988 
(1963); R.Z. Sagdeev and A.A. Galeev, Nonlinear Plasma 
Theory (Benjamin, New York, 1969); M.L. Goldstein, As- 
trophys. J. 219, 700 (1978). 

[4] F. Malara, L. Primavera and P. Veltri, Phys. Plasmas 7, 
2866, (2000). 

[5] R.H. Cohen and R.M. Kulsrud, Phys. Fluids 17, 2215 
(1974); V.M. Nakariakov, L. Ofman and T.D. Ar- 
ber, Astron. Astrophys. 353, 741 (2000); E. Verwichte, 
V.M. Nakariakov and A. Longbottom, J. Plasma Phys. 
62, 219, (1999). 

[6] J. Hayvaerts and E.R. Priest, Astron. Astrophys. 117, 
220 (1983). 

[7] F. Malara, L. Primavera and P. Veltri, Astrophys. J. 459, 
347, (1996). 

[8] L. Ofman and J.M. Davila, J. Geophys. Res. 100, 



23413, (1995); Astrophys. J. 476, 357, (1997); S. Poedts, 
G. Toth, A.J.C. Belien and J.P. Goedbloed, Solar Phys. 
172, 45, (1997); R. Grappin, J. Leorat and A. Buttighof- 
fer, Astron. Astrophys. 362, 342 (2000); I. De Moortel, 
A.W. Hood and T.D. Arber, Astron. Astrophys. 354, 334 
(2000). 

[9] V.M. Nakariakov, B. Roberts and K. Murawski, Solar 
Phys., 175, 93 (1997); Astron. Astrophys. 332, 795 
(1998). 

[10] I.R. Mann, A.N. Wright and P.S. Gaily, J. Geophys. Res. 

100, 19441 (1995); M. Goossens and A. De Groof, Plas. 

Physics 8, 2371 (2001) 
[11] G.J.J. Botha, T.D. Arber, V.M. Nakariakov and 

F.P. Keenan, Astron. Astrophys., 363, 1186 (2000); 

D. Tsiklauri, T.D. Arber and V.M. Nakariakov, Astron. 

Astrophys., in press (2001) 
[12] T.D. Arber, A.W. Longbottom, C.L. Gerrard and 

A.M. Milne, J. Comput. Phys., 171, 151 (2001). 
[13] M. Kaviany, in The Handbook of Fluid Dynamics (Ed. by 

R.W. Johnson) (CRC Press, Boca Raton: Fla., 1998), p. 

21-34. 



